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ABSTRACT 

The dynamics of the wind- wind collision in massive stellar binaries is investigated using 
three-dimensional hydrodynamical models which incorporate gravity, the driving of the 
winds, the orbital motion of the stars, and radiative cooling of the shocked plasma. 
In this first paper we restrict our study to main-sequence O+O binaries. The nature 
of the wind-wind collision region is highly dependent on the degree of cooling of the 
shocked plasma, and the ratio of the flow timescale of the shocked plasma to the orbital 
timescale. The pre-shock wind speeds are lower in close systems as the winds collide 
prior to their acceleration to terminal speeds. Radiative inhibition may also reduce the 
pre-shock wind speeds. Together, these effects can lead to rapid cooling of the post- 
shock gas. Radiative inhibition is less important in wider systems, where the winds 
are accelerated to higher speeds before they collide, and the resulting collision region 
can be largely adiabatic. In systems with eccentric orbits, cold gas formed during 
periastron passage can persist even at apastron, before being ablated and mixed into 
its surroundings and/or accelerated out of the system. 

Key words: hydrodynamics - shock waves - stars: binaries: general - stars:early-type 
- stars: mass loss - stars: winds, outflows 



1 INTRODUCTION 

The collision of the winds in massive star binary sys- 
tems can produce a rich array of phenomena. Energetic 
particles are a ccelerated at the wind-wind collision region 
(WCR -see eg.lEichler fc Usovll99H;lDougherty et al.ll2003l ; 
IPittard et all 120061 ; IPittard fc Dougherty! l200rj ), producing 
non -thermal radio emission via the synchrotron process (see 
e.g. iDoughertv et al]|2005l ). and X-ray and 7-ray emission 
via the inverse Compton and pion decay processes which 
should extend to T eV energies (|Pittard fc Dougherty! I2OO6I: 
iReimer et al.ll2006h . Dust can be formed either episodic ally 
producing transient infra-red outbursts (|Williamslll99r3 ), or 
conti nuously producing s piral-shaped structures on the sky 
(e.g. iTuthill et aT] Il999l : (Marchenko et al.l 120021 ). Emission 
line profiles from the WCR show variable doppler shifts and 
broadening as sight lines relative to the WC R change (e.g. 
lLuhrslll997i ; lHenlev. Stevens fc Pittardll2003h . Shock heating 
of the winds produces bright thermal X-ray emission, which 
may e xhibit signs of non- equilibrium ionization in wide sys- 
tems |Pollock et al.ll2005t ). 

The hydrodynamics of colliding wind binaries (CWBs) 
have usually been investigated using two-dimensional 
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(axis ymmetric) models with termin a l velocity winds 
(e.g. IStevens. Blondin fc Pollock! Il992l : IPittard fc Stevens! 
Il997l ). 3D simulatio ns whic h inclu d e orbital effect s have 
been presented by IWalderi (| 19981 ). IPittard! i| 19991 ). and 
iLemaster et al.l (|2007 ). In close binaries, the winds col- 
lide before they reach terminal velocity, and their dynam- 
ics may be altered by the companion star's radiation field 
JSteyens fc Pollock! Il99i IGavlev. Owocki fc Cranmerlll997l . 
Il999l ). The driving of the winds and the dynamical effect of 
the companion's gravitational and radi ation fields has bee n 
included in 2D models of V444 Cy gni (IGavlev et al.lll997l ). 
l Orionis dPittardl 19981). n Carinae (|Pittard et al.lll99Sl ), and 
Sanduleak 1 (|St-Louis et al.ll2005l ). The only publi shed 3D 
simu lation with radiative driving is of V444 Cygni (|Pittardl 
Il999l ). Other approaches based on solutio ns to the ram- 
pressu re b alance have been presen ted by lAntokhin et al.l 
|2004l ) and lParkin fc Pittardl feoOSl ). 

In this work we examine the hydrodynamic properties 
of the wind-wind collision in short period O+O binaries. 
We perform 3D calculations, so that orbital effects can be 
included, and calculate the radiative driving of the winds, so 
that the winds collide at realistic velocities and some effects 
of the companion's radiation field can be accounted for. In 
some of the simulations the winds collide at only a fraction 
of their terminal velocities, and the lower postshock temper- 
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atures and higher postshock densities which result enhance 
the radiative cooling of the postshock gas. We examine the 
dynamics of the wind- wind collision in systems with circular 
and eccentric orbits, and with equal and unequal winds. In 
Section [2] we describe the details of our simulations and in 
Section Owe present our results. Section U summarizes and 
concludes this work. 



2 THE NUMERICAL SETUP 

2.1 Details of the hydrodynamics 

Our three-dimensional simulations were conducted on a 
cartesian grid using Eulerian hydrodynamics with piecewise 
parabolic interpolation of the fluid variables. The code solves 
a Riemann problem at each zone interface to determine the 
time-averaged values at the zone faces, and then solves the 
equations of hydrodynamics: 

dp 



Table 1. Assumed binary parameters for the models investigated. 
The semi-major axis is 34.26 Rq in model cwbl, 76.3 Hq in 
models cwb2 and cwb3, and 55 R© in model cwb4. 
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where e = u 2 /2+e/p is the total specific energy, p is the mass 
density, e is the internal energy density, P is the pressure, 
and T is the temperature. We adopt an ideal gas equation 
of state, e = Pj (7 — 1), with 7 = 5/3 as the ratio of specific 
heats, f is the force per unit mass and includes gravity and 
radiative driving terms. 

The radiative cooling term, A(T), is calculated from the 
MEK AL thermal plasma code (|Mewe et all Il995l : iKaastral 
11992) 1 distributed in XSPEC (vll.2.0). The plasma is as- 
sumed to be in collisional ionization equilibrium (but see 
Sec. I2.2|l . The temperature of the pre-shock stellar winds 
is assumed to be maintained at ~ 10 4 K through pho- 
toionization heating by the stars. Gas in the WCR which 
rapidly cools is prevented from cooling below this tempera- 
ture. Because of the high Mach numbers involved, the den- 
sity contrast of the hot plasma and cooled regions can be 
very high (the density contrast across an isothermal shock 
of Mach number M is 7A/ 2 ). In reality, magnetic pressure 
may halt the compr ession before this ratio is reached (e.g. 
iKashi fc S okcr 200?]), but this is an additional complication 
which is not considered here. The code also contains several 
advected scalars which allow tracking of which wind mate- 
rial is in which cell, and the ionization age and temperature 
equilibration of the postshock gas (see Section f2.2p . 

The body forces acting on each hydrodynamic cell are 
the vector summation of gravitational forces from each star, 
and continuum and line driving forces from the stellar ra- 
diation fields. The computatio n of the line acceleration is 
based on a local ISobolevI (|l960l) treatment of the line trans- 
port, following the standard Castor. Abbott fc Kleinl (1 19751 . 
hereafter CAK) formalism developed for single OB winds, 
au gmented by the finite disk c o rrecti on factor developed 
by IPauldrach. Puis fc Kudritzkil (|l986l ). The vector radia- 
tive force per unit mass, g rad , is computed from an integral 
of the intensity /(n) times the projected velocity gradient 
along directional vectors fi within the solid angle covering 
the stellar disk, 



Model 


Stars 


Period 


Eccentricity 


Wind mtm. 






(d) 


(e) 


ratio (rj) 


cwbl 


O6V+O6V 


3 


0.0 


1 


cwb2 


O6V+O6V 


10 


0.0 


1 


cwb3 


O6V+O8V 


10.74 


0.0 


0.1 


cwb4 


O6V+O6V 


6.1 


0.36 


1 



Table 2. Assumed stellar parameters for the models investigated. 
Note that the value of Q is fairly constant, and shows less varia- 
tion than k in the two stellar models constructed for this work. 



Parameter /Star 


O6V 


O8V 


Mass (M Q ) 


30 


22 


Radius ( R© ) 


10 


8.5 


Effective temperature (K) 


38000 


34000 


Mass-loss rate ( Mq yr -1 ) 


2 x 10~ 7 


10- 7 


Terminal wind speed (kms - 1 ) 


2500 


2000 


CAK a 


0.57 


0.52 


CAK k 


0.12 


0.24 


Q 


258 


308 


C J pVth 


a 

hdfl. 



(4) 



The integration over solid angle is performed with 8 direc- 
tional vectors, a and k are the standard CAK parameters, a e 
is the specific electron opacity due to Thomson scattering, 
and v t h is a fiducial thermal velocity calculated for hydro- 
gerQ. Shadowing by the companion star is accounted for in 
our calculations, and in such cases only the visible part of 
the stellar disk contributes to the radiative driving force. 
The line driving is set to zero in cells with temperatures 
above 10 6 K, since this plasma is mostly ionized. Further 
de tails concerning the line force calculations can be found 
in lGavlev etHI (| 19971 ). 

A companion star influences the driving of a wind in a 
number of distinct ways. Very close to the surface of the star, 
reflection of the companion's starlight from the stellar photo- 
sphere can actually enhance the wind driving and mass loss, 
and the angular asymmetry in the line-scattering probabil- 
ity can induce significant tilting of t he direction vector o f the 
wind away from the surface normal (|Gavlev et al II 19991 ) . Re- 
flection effects become less important further from the star. 
At this point the main effect of the companion's radiation 
field is to reduce the net radiative flux, slowing the remaining 
acceleration and producing a lower terminal velocity - this 



1 While we use the standard CAK a and k description, we note 
that the k parameter does not represent a physically meaningful 
quantity in its own right, and has an artifical dependence on v^. 
A more mea ningful descrip tion of the strength of line driving 
was given by iGavlevI l ll995l ). who introduced the Q parameter. 
a and k values can easily be converted into Q values, and are 
noted alongside the a and k values for each of our model stars in 
Tabled 
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Table 3. Details of the hydro-dynamical grid of each model. The 
length of the sides of the grid in model cwb4 corresponds to 3.2 — 
6.86 D Bep as the stars progress in their orbit from apastron to 
periastron. 



Model 


No. of cells 


Length of sides 


Resolution 






(R ) (a) 


(R©) 


cwbl 


480 3 


240 7 


0.5 


cwb2 


456 3 


570 7.46 


1.25 


cwb3 


456 3 


570 7.46 


1.25 


cwb4 


480 3 


240 4.34 


0.5 



effect is known as radiative inhibition (|Stevens fc Pollockl 
1 19941 ). The gravity of the companion star also has an effect - 
it alters the net force on the wind, and can distort the shape 
of the other star, which may lead to gravity darkening of the 
stellar surface and thus again affect the wind. 

Outflow boundary conditions are used for all of the sim- 
ulations. The fluid values on the grid are initialized by map- 
ping ID solutions of the CAK equations onto the grid. The 
wind velocities are modified to account for the velocity of 
each star on the grid (i.e. v = v w (r) + We assume that 
the stars remain closely spherical and are not subject to 
significant deformations due to the companion's gravity. To 
generate the winds we re-initialize density, pressure and ve- 
locity values within a shell of 3-cells thickness around each 
star at the start of every timestep. This procedure means 
that the stars maintain a constant mass-loss rate, and have 
a wind which is initially directed radially outwards. Because 
our models do not resolve the wind acceleration regions very 
close to the stars, we are unable to investigate reflection ef- 
fects on the mass-loss and wind dynamics. Instead the focus 
is on exploring the latter stages of the wind acceleration and 
the effects this has on the resulting wind-wind collision. 

Our treatment of the wind initiation also means that 
the winds effectively feel no radiative inhibition effects un- 
til outside the 3-cell remap radius. While this is not en- 
tirely satisfactory, our attempts to self-consistently acceler- 
ate the winds from the stellar surfaces have to date always 
resulted in unacceptable wind structure due to a "staircase" 
effect introduced by the grid zones which straddle the sur- 
faces of the stars. This problem is very much reduced in 2D 
r — 6 simulations of single stars when the sonic point of the 
winds is resolved (see e.g. lOwocki et al.|[l994l ). To resolve 
this problem in 3D simulations with more than one star 
may require adaptive-mesh-refinement (AMR) calculations, 
and/or sophisticated grid geometries. For the time being we 
proceed with the above limitations, with the aim of reduc- 
ing/eliminating this problem in future work. 

We specify the initial boundary between the winds as a 
flat 2D plane which is normal to the line of centers through 
the stars and passes through the stagnation point where 
momentum balance is achieved. This position is calculated 
assuming the winds instantaneously reach terminal velocity. 
The simulations are evolved until the effects of the initial 
conditions have been driven off the grid, following which we 
start our analysis. For numerical reasons the orbital plane 
was aligned with the yz plane on the grid, with the x-axis 
perpendicular to this. 



2.2 Postshock equilibration 

There is now strong evidence from studies of supernova 
remnants that the postshock thermalization of electrons 
lags behind that of ions in high speed collisionless shocks, 
with the ratio of postshock electron to proton tempera- 
ture, T /Tp, su bstantially below unity for shock velocitie s 

> 1000kms _1 (|Rakowskill2005l ; I van Adelsberg et al.|[2008l ). 
Temperature equilibration through Coulomb collisions and 
plasma instabilities then occurs some distance downstream. 
One may expect these findings to be of relevance to collision- 
less shocks in general and thus also to CWBs. In fact, the 
postshock thermalization of electrons has already been ex- 
am ined in the prototypical colliding winds binary WR 14 
by IZhekov fc Skinner! (|2000l ) and by IPollock et all (|2005l) . 
and can naturally account for t he softer-than-expec ted X- 
ray continuum emission (see also lPittard et aTll2006l ) . 

Despite an increased understanding of the need to 
consider slow electron heating, the physics of the heat- 
ing process in the shock layer remains poorly understood. 
For instance, the appropriate value of T e /T p immediately 
downstream of the subshock is not clear, with observations 
yielding a minimum value of ~ 0.03 at shock velocities 

> 1500 km s -1 which is greater than the theoretical mini- 
mum (me/jrip) by two orders of magnitude, but smaller than 
the predictions of some coll isionless heating models (e.g. 
ICargill fc Papadopoulosll 19881 ). Faced with this situation, we 
have therefore assumed in our models that T e /T — 0.2 im- 
mediately post shock, where T is the m ean plasma temper- 
ature (see also lZhekov fc Skinner] [2OO0I ), The chang e in the 
down stream electron temperature is then given by (|Spitzerl 
[1973) 



d(T c /T) 
dt 



T e \ In A / T c 



= 3.8xl0- 12 (n c +n i )(l-^)^ 



VlO 7 K 



-3/2 



where the numerical constant assumes cosmic abundance s 
and In A « 30 (cf. iBorkowski. Sarazin fc Blondinl [19941 ) . 
Equipartition requires that the product of the postshock 
timescale and the density exceed nt > 3 x 10 11 cm -3 s. 
Clumping in the win ds will also a ffect the electron ther- 
malization timescale (|Pittard|[2007l) . Eq. \5\ is solved along 
with the other hydrodynamical equations. 

The timescale for the postshock ionization to approach 
equilibrium will also be an important factor in some CWBs. 
For instance, non-equilibrium ionization may explain why 
line-p rofile models which assume rapid ionization equilib- 
rium l|Henlev et all 120031 ) are unable to reproduce the ob- 
served correlation of X-ray line widths with i onization po- 
tenti al in the Wolf-Rayet system 7-Velorum (I Henley et alj 
2005j). Although the postshock ionization depends on the 
thermal history of the plasma, to zeroeth order the ion- 
ization is independent of the tempe rature histor y and the 
specific ionic species in the plasma (|Masai|[l994l ). and can 
be characterized by the so-called ionization age, n e t. For 
example, the ionization of oxygen and iron in plasma at 
T = 10 6 5 K as a func t ion of n e t is shown in Figs. 1 and 2 
of iHughes fc Helfandl (| 19851 ). The oxygen reaches ioniza- 
tion equilibrium at n c t ~ 10 12 ' 5 cm -3 s, while the iron 
reaches equilibrium at n e t 10 11 ' 5 cm~ 3 s. At higher tem- 
peratures, iron is ionized to higher stages, and requires a 
greater ionization age to reach equilibrium. However, in gen- 
eral a plasma will be close to ionization equilibrium when 
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Table 4. Some key parameters of each model. In model cwb3 two values are given for many of the parameters - the first corresponds 
to the primary (06V) star/wind, and the second to the secondary (08V) star/wind. For model cwb4 values at periastron and apastron 
are noted. v OI ^ is the orbital speed of the stars and u w is the preshock wind speed along the line-of-centres, both in kms -1 . i> w /^s is 
the ratio of the preshock wind speed along the line-of-centres to the wind speed at the same radial distance from a single-star solution. 
Smaller values indicate more effective radiative inhibition. The ratio v^^/v-w affects the aberration angle, a j,, of the WCR. The values 
quoted for ab are measured directly from the simulations (not calculated from the simple formula in Scction l2.3t . and noted in degrees. 
The degree of downstream curvature of the WCR in the orbital plane is specified by a COT , where the curvature is assumed to trace an 
Archimedean spiral which in polar coordinates is described by r = a cav Q. The value of a C or corresponds to the downstream distance (in 
units of d Bcp ) along the WCR for each radian of arc it sweeps out in the orbital plane. Smaller values indicate tighter curvature. The 
leading and trailing arms of the WCR in model cwb3 display differing degrees of curvature, so the value quoted for this model is a rough 
average. For model cwb4 the quoted values roughly represent the "instantaneous" curvature at periastron and apastron. \ is the ratio of 
the postshock cooling to flow timescales, and is calculated using the pre-shock wind speed, i> w , as measured from the simulations. \ < 1 
indicates that the shocked gas rapidly cools, while x ^ 1 indicates that the plasma in the WCR remains hot as it flows out of the system. 



Model v orh v„ Vw/vs v oy . h /v w ab a COI \ 



cwbl 


290 


730 


0.58 


0.40 


17 


3.5 


0.34 


cwb2 


225 


1630 


0.90 


0.14 


3-4 


6.5 


19 


cwb3 


152,208 


1800,1270 


0.90,0.88 


0.08,0.16 


2 


4.5 


28,14 


cwb4 


334-156 


710-1665 


0.57-0.91 


0.47-0.09 


21-4 


3-10 


0.34-19 



n c t > 10 12 ' 5 cm -3 s. We therefore advect a scalar variable 
in the hydrodynamical code to track the postshock value of 
n c t. 



2.3 Models investigated 

In this work we do not attempt to model particular sys- 
tems. Instead the aim is to achieve an understanding of 
how the dynamics of the collision region depend on some 
key parameters. We have therefore computed 4 different 
models with a range of binary and stellar parameters, as 
summarized in Tables [T] and [2] In the present work we 
consider systems with dwarf O-stars, so that tidal distor- 
tions are minimized, and in this first study can be ig- 
nored. The stellar para meters are consistent with recen t 
observations presented inlRepolust. Puis fc Herrerol d2004h. 



iMartins. Schaerer fc Hillierl ( 20051 ). iBouret, Lanz fc Hillierl 
(|2005h and lFullerton et al. M2006J ). in which clumping within 
the winds is considered. The CAK parameters (a and k) 
which are needed to drive the desired winds from the stars 
are noted in Table along with the value of Q calculated 
from them. We assume that the mass-loss from the stars is 
isotropic. Solar abundances are assumed for each star/wind, 
and all simulations are performed in the centre of mass 
frame. 

The nature of the WCR is largely governed by the ra- 
tio of the cooling timescale for the shocked gas to the dy- 
namical timescale for it to flow out of the system, which 
is given by x — Vgdi2/M-7, where vg is the preshock ve- 
locity in units of 1000 kms -1 , di2 is the stellar separa- 
tion in units of 10 12 cm, and M -7 is the stellar ma ss-loss 
rate in units of 10" 7 M yr" 1 (|Stevens et al.l 1 19921 1. An- 
other key parameter governing the nature of the WCR is 
the ratio of the orbital speed to the pre-shock wind speed, 
^orb/«w- This ratio affects the aberration angle of the WCR 
[8 ~ tan -1 (-!;orb/«w)] and its degree of downstream curva- 
ture due to coriolis forces. 

In this work we consider systems with identical or 
similar wind mome nta. This means that radiative braking 
(|Gavlev et alj|l997h effects are negligble, and in all cases a 



normal wind-wind pressure balance exists. In future work we 
will extend our investigation to systems with higher wind 
momentum ratios, to study the effects of braking and the 
dynamics of the WCR when the more powerful wind over- 
whelms the weaker wind and the WCR collapses onto the 
surface of one of the stars. 

In model cwbl we suppose that two identical 06V stars 
move around each other in a circular orbit with a period of 3 
days. The stellar separation is 34.26 Rq, and each star has 
an orbital velocity v OT \, = 290 kms -1 . A rough estimate of 
X can be obtained by determining the speed of the winds at 
the stagnation point where momentum balance is achieved. 
Ignoring for the moment effects due to the orbital motion of 
the stars and the presence of two radiation fields, one finds 
that the winds accelerate to a speed of 1250 kms -1 prior 
to their collision at the stagnation point. This implies that 
X = 2.9, and the postshock gas would not be expected to 
cool back to T ~ 10 4 K until it had left the central regions 
of the system. In fact, we shall see that once the simula- 
tion starts, the reduction in the net radiative flux due to 
the presence of the companion star reduces the accelera- 
tion of eac h wind towards its comp anion. This "radiative 
inhibition" (|Stevens fc Po llock 199J) reduces the pre-shock 
velocities below 1250 kms -1 , which increases the cooling 
in this model (see Table 2] which summarizes some of the 
salient features directly obtained from the hydrodynamcial 
simulations). Neglecting radiative inhibition again, we also 
find that v OT b/v w = 0.23, indicating that the WCR will dis- 
play a reasonably large aberration angle and downstream 
curvature due to coriolis forces. These effects become more 
pronounced if there is a significant reduction in wind speeds 
along the line-of-centres due to raditive inhibition. Note also 
that unphysical masse s are required in some of the models of 
iLemaster et al.l (|2007t ). including their main large-box sim- 
ulation C2.5. 

Model cwbl is s imilar to many r eal systems, in- 
cluding DH Cep (see Linder et al.l 20071. and reference s 



therein), HP 16 5052 (lArias et al 



and HD 159176 (|De Becker et al.ll2004bl : 



[200% 



Lin der et al]|2007l) , 



Linder et afll2007h . 



All of these systems have near identical main-sequence stars 
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Figure 1. Density plots of models cwbl (a), cwb2 (b), and cwb3 (c). The colour scale is logarithmic, spanning 10~ 19 gcm~ 3 (blue) to 
10 -11 gem -3 (red). Panel (a) has sides of length 240 Rq, while panels (b) and (c) have sides of length 570 Rq. T he black, red and 
green lines marked on the figures show the position of the contact discontinuity as calculated by the model described in lParkin fc Pittardl 
l|2008h - see the text for further details. 




Figure 2. a) 3D density contours from model cwbl. The orange contour shows the surface at a temperature of 3 X 10 K, while the red 
contour shows the surface at a density of 10~ 14 gcm~ 3 . The density contour highlights the stellar positions, plus the densest regions 
of the WCR (the other star is hidden behind this). The WCR readily breaks up into small dense clumps on the trailing edge of each 
arm. These clumps are often surrounded by high temperature bow shocks as revealed by the temperature surface, b) As Fig. I^a) but 
for model cwb2. The shocks trace a double-helix, in a similar manner to the DNA molecule. 




Figure 3. Velocity plots of models cwbl (a), cwb2 (b), and cwb3 (c) in the orbital plane. The colour scale is linear, spanning kms 
(blue) to 4500 kms -1 (red). 
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Figure 4. Temperature plots of models cwbl (a), cwb2 (b), and cwb3 (c) in the orbital plane. The colour scale is logarithmic, spanning 
^ 10 5 K (white) to 10 s K (pink). 



of spectral type 06-07, and circular or near-circular orbits 
with periods near 3 days. Hence the hydrodynamics of and 
emission from the WCR in model cwbl will be a reasonable 
approximation to the situations in these systems. In some 
of these systems it may be difficult for the winds to attain 
a normal ram-ram pressure balance, and the stronger wind 
may push the WCR onto the surface of the star with the 
weaker wind. Detailed simulations of each individual sys- 
tem are needed to ascertain whether, in fact, this occurs. 
Reflection and distortion effects may also well be important 
in these systems, but are beyond the scope of the present 
work. 

Model cwb2 uses the same stars as model cwbl, but 
increases the orbital period to 10 days (the stellar separa- 
tion becomes 76.3 R©). The pre-shock collision speed on the 
line of centers (neglecting radiative inhibition) increases to 
1970 kms -1 , as the winds now have more room to accel- 
erate. The higher collision speed leads to higher postshock 
temperatures, and together with the increased stellar sepa- 
ration produces lower postshock densities. The result is that 
cooling within the WCR is dramatically reduced {% m 40), 
and the postshock gas behaves largely adiabatically. The re- 
duced orbital speed and higher pre-shock wind speeds also 
means that the aberration and curvature of the WCR will 
be smaller in this model than in model cwbl. Taking ac- 
count of radiative inhibition we again find slightly increased 
cooling and orbital effects (see Table [3J|. 

Model cwb2 is similar to HD 93161A, an 08V + 09V 
syste m with a circular orbit and an orbital period of 8.566 
days |Naze et alj|2005h . albeit with slightly more massive 
stars and powerful winds. Another sy stem with not too 
dissim ilar properties is Plaskett's star (|Linder et al.|[200?3 . 
120081 ). though this object contains stars which have evolved 
off the main sequence. The wider separation of the stars in 
these systems means that the dominant dynamical effect of 
the companion starlight will be radiative inhibition, with 
reflection effects minimized. 

In model cwb3, we examine the interaction of unequal 
winds in a hypothetical O6V + 08V binary. We keep the 
same stellar separation as model cwb2, which gives a pe- 
riod of 10.74 days from the total system mass of 52 Mq. 
The orbit is circular, with the O8V star further from the 
center of mass. An initial estimate of the preshock wind 
speeds can be obtained by considering the location of the 



stagnation point. The terminal speed wind momentum ra- 
tio, r\ — M2«oo,2/Afi« 0Oj i = 0.4, indicating that the stagna- 
tion point occurs at a distance of 0.39 a from the O8V star. 
At this distance the secondary wind speed is 1570 km s -1 , 
giving \2 ~ 32. The primary wind collides at higher speed 
(2070 kms^ 1 ), and its postshock plasma is slightly more adi- 
abatic (xi ~ 50). Again, radiative inhibition reduces the 
preshock wind speeds and postshock cooling parameters be- 
low these simple estimates (see Table |3J. 

Model cwb4 explores the effect of an eccentric or- 
bit which takes the stars through separations of 34.26 — 
76.3 R0 (i.e. the separations of the stars in the circu- 
lar orbits of models cwbl and cwb2). The required ec- 
centricity is e = 0.36, and the resulting orbital period is 
6.1 days. This model allows us to investigate whether the 
WCR properties at a specific stellar separation are similar 
to those obtained with a circular orbit. Some well-known 
O+O binaries with eccentric orbits include (in order of in- 
creas ing orbital period) HP 152248 (e = 0.127: ISana et all 
|2004 ). HP 9 3205 (e = 0.46: iMorrell et alJkoOlh , HP 93403 



(e = 0.234: iRauw et all I2002D. CygOB2#8A (e = 0-24 ; 
|Pe Becker. Rauw fe Manfroidl 12004: |Pe Becker et al.l [iooeh 



and f-Orionis fe = 0.764: iBagnuolo et alj|200ll ). 

The hydrodynamical grid is half-cubic in all simula- 
tions, and is reflected in the orbital plane for further anal- 
ysis. The grid is large enough to capture, for example, the 
majority of the X-ray emission from each model. Petails of 
the grids are noted in Table [3] As a compa rison, most of the 
simulations run bv lLemaster et all (|2007i ) used a grid with 
sides of length 2.5 a, with their large grid simulations having 
sides of length 6.25 a. 



3 RESULTS 
3.1 Model cwbl 

Fig. QJa) shows a density snapshot in the orbital plane of 
model cwbl. The postshock gas in the WCR between the 
stars has cooled into a thin, dense sheet, which is sub- 
ject to the non-linear t hin-shell instability ljVishniadll99i : 
iBlondin fe Markslll995h . Clearly our initial estimate of \ in 
Section [2.31 was too high. Examination of the dynamics re- 
veals that the winds in fact collide at « 730 km s -1 , rather 
than the 1250 km s" 1 obtained from considering the accel- 
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oration of the wind of an isolated star. This is caused by a 
significant reduction in the acceleration of the stellar winds 
towards the centre of mass of the system due to strong radia- 
tive inhibition (in reality reflection effects could also be im- 
portant, though this process is not examined in this work). 
The cooling parameter, x> is reduced to ~ 0.3, which gives 
rise to the severe radiative cooling of the postshock gas seen 
in Fig. [Ha). 

As previously noted, the nature of the WCR is also de- 
termined by the ratio of the pre-shock wind speed to the or- 
bital speed of the stars, since this affects its angle of aberra- 
tion and downstream curvature. The aberration angle of the 
WCR in Fig. QJa) is ~ 17°, which compares to an expected 
value of 22° if v m = 730 km s -1 is used in the formula in 
Section \l. 31 That the expected value is overestimated is en- 
tirely anticipated, of course, since the pre-shock wind speed 
is higher off-axis than it is along the line-of-centres, so that 
using speeds along the line-of-centres biases the calculation 
to larger angles. 

Downstream the WCR spirals around the stars as a re- 
sult of their orbital motion. Analyses of the so-called "pin- 
wheel nebulae" have compared such morphologies to an 
Archimedean spi ral, which in polar coordinates is specified 
byr 9 Csee lTuthill et alj|2008l . and references therein) . 

The curvature of the spiral is governed by the distance which 
the plasma in the WCR moves during one orbital period. 
The gas in the denser parts of the WCR near the right 
edge of the grid shown in Fig. [TJa) has a speed of about 
2050 kms -1 . Hence in 3 days this travels a distance of 
3.55 au (763 Rq). The minimum distance from the cen- 
tre of the grid to its edge is 120 Rq , so we expect the WCR 
to curve through 16 per cent of a full rotation, or about 57°. 
Fig. [TJa) actually displays greater curvature than this esti- 
mate, though this is expected because the speed of gas near 
the apex of the WCR is much slower than further down- 
stream. 

Overplotted on Fig. [Ha) is the contact discontinuity 
in the orbit al plane as calcula t ed by the dynamical model 
described in iParkin fc Pittardl (|2008l ). The model assumes 
that the winds collide at terminal velocity, which is ob- 
viously not the case here. Nevertheless, the model is able 
to match remarkably well the position of dense gas in the 
WCR, if the wind speeds are assumed to be 2000 kms - 
and the flow in the WCR behaves ballistically once acceler- 
ate d to 75 per cent of its te rminal speed (for further details 
see IParkin fc Pittardl 120081 ). Adopting a 2500 kms -1 ter- 
minal wind speed in the dynamical model results in slightly 
less aberration and downstream curvature, while specifying 
a 750 kms -1 terminal wind speed results in much greater 
aberration and a much tighter downstream curvature com- 
pared to the hydrodynamical models. Clearly, there is scope 
for the dynamical model to be "calibrated" against full hy- 
drodynamical models. A major advantage of the dynamical 
model over full hydrodynamical calculations is its low com- 
putational cost. 

From movies of the simulation it is clear that the radia- 
tive overstability occurs in model cwbl. High Mach num- 
ber shocks are susceptible to a global oscillation if the 
slope, a, of the cooling function (where the cooling func- 
tion, A(T), is approximated as Apr"), is less than a criti- 
cal value (i.e. a < 0.4 - see e.g. IStrickland fc Blondinlll995l ; 
iPittard et ai]|2005f ). In this model the postshock gas tem- 



perature T ~ 7 x 10 K corresponds to a local slope a ~ 
on the cooling function. This is less than the critical value 
a cr = 0.4, below which the overstability occurs. As the gas 
cools below ^3xl0 6 K, very low values of a are encountered 
(see IPittard et all 120051 ) . Mass- loading (from clump abla- 
tion and evapor ation) can moderate or even prevent ther- 
mal instabilities (|Pittard. Hartquist fc Ashmorel2003bh , but 
clearly is not a significant effect here. 

Further downstream, the leading edge of the WCR be- 
comes shadowed by upstream regions of the WCR, lead- 
ing to very low post-shock densities, and very weak oblique 
shocks. In contrast, the trailing edge of the WCR is very 
much denser, and remains much more normal to the wind. 
The WCR remains denser on the trailing side because of 
the high inertia of the dense gas, which causes it to respond 
slowly to the movement of the stars in their orbits. The dense 
sheet of cooled material between the stars is broken up into 
numerous small clumps in the downstream flow on the trail- 
ing side of the WCR, as shown in Fig. [2ja) which shows 3D 
isosurfaces of the density and temperature in model cwbl. 
In contrast, the trailing edge of the WCR is far smoother 
and less structured. 

Fig. [3ja) shows the wind speeds in the orbital plane. 
The driving of the winds off the back of the stars produces 
wind speeds which are lower in the centre of mass frame 
than off the front of the stars. Inspection of the model also 
reveals that the winds accelerate more slowly when in the 
shadow of one of the stars. Remarkably, the maximum speed 
of material on the computational grid in the orbital plane is 
~ 4650 kms - , almost double the terminal velocity of the 
single star wind. However, these speeds occur in extremely 
rarefied gas within the WCR and are atypical of the highest 
speeds seen. It is unlikely that such high speeds would be 
easily detected since they are tied to such small amounts 
of material. A more representative speed of the high veloc- 
ity gas is 3000 — 3500 kms -1 (see Fig. (7]d) - such speeds 
are typically seen in the shadow of the leading shocks of the 
WCR. The direction of the lowest wind speeds is towards the 
top left and bottom right of the density plot in Fig. [TJa), 
where the typical wind speed is ~ 2750 km s -1 . The winds 
are slower in these directions because of the oppositely di- 
rected stellar motions. The winds are still accelerating in 
these regions, albeit gradually. 

Fig. |4ja) shows the temperature in the orbital plane of 
model cwbl. The expected postshock temperature behind a 
high Mach number shock which is normal to the preshock 
flow is given by T = jg^v 2 K, where m is the average 
mass per particle. Given the observed preshock wind speeds 
along the line of centres we expect a postshock temperature 
T = 7.2 x 10 6 K, which is indeed what is seen. 

In wide CWBs where the winds collide at their terminal 
speeds the maximum postshock temperature is produced at 
the apex of the WCR where the shocks are normal to the 
flow. For a given preshock speed, oblique shocks are less effi- 
cient at thermalizing the flow and produce lower postshock 
temperatures. However, in systems where the winds collide 
prior to reaching their terminal speeds the hottest tempera- 
tures in the WCR can in fact be obtained far downstream of 
the WCR apex, as we now explain. The increased distance 
from the stars to the shocks means that the winds have more 
time to accelerate and collide at higher speed. While these 
increased speeds may be offset by the increasing obliquity of 
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Figure 5. Density (top), temperature (middle), and speed (bottom) slices of model cwbl in the orbital plane (left), and below the 
orbital plane (middle and right). The density and temperature colour scales are the same as in Figs. [T] and [4] while the colour scale for 
the speed plots is the same as the temperature scale though is linear and covers the range — 4500 km s~ 1 (i.e. stationary/low speed gas 
is white/blue, while the high speed gas is red). 



the shocks, in cases where the WCR is strongly radiative and 
thus susceptible to kinking by thin shell instabilities, some 
parts of the WCR surface (and thus the shocks) can be nor- 
mal to the preshock flow. This combination of high preshock 
speeds and normal shocks can result in the hottest plasma 
temperatures actually occuring far downstream. This is in 
fact the case in model cwbl, where temperatures reach as 
high as 4.5 x 10 7 K at the apex of bowshocks around down- 
stream clumps, indicating a relative speed of ~ 1800 kins -1 
between the dense clumps and the preshock winds at such 
points (though the typical temperature in these bowshocks 
is somewhat lower - see Fig. [7p). Note that the maximum 
temperature is actually higher than the corresponding value 
from model cwb2 (see Section f3 . 2 f> . despite the latter mod- 
els wider stellar separation and higher preshock wind speeds 
along the line-of-centres (at first glance Fig.JTfb) would seem 



to contradict this statement, but this is only because the or- 
dinate as plotted does not extend to low enough values). 

In model cwbl, the densities between the stars are suffi- 
ciently high for the shocks to be collisional. Hence, the elec- 
tron and ion temperatures rapidly come into equilibrium. 
The post-shock density (prior to cooling) along the line-of- 
centres is n ~ 10 10 cm -3 (this increases to a maximum 
of 2.4 x 10 12 cm~ 3 after cooling - see Fig. [7^). Hence the 
electron and ion temperatures reach equilibrium about 30 
seconds after passing through the shocks. The equilibration 
timescale is longer for gas shocked downstream of the WCR 
apex, but still very rapid. 

Similarly, the ionization age of the shocked plasma, n e t, 
exceeds 10 13 cm -3 s near the apex of the WCR, indicating 
that the post-shock gas is in collisional ionization equilib- 
rium. Lower ionization ages are seen further downstream 
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Figure 6. Density (left), temperature (middle), and speed (right) slices of model cwbl in the plane orthogonal to the orbital plane 
containing the stars (top), and in the plane orthogonal to this and the orbital plane (bottom). The colour scales are the same as in Fig. [5] 



(e.g. in the bowshocks surrounding the protrusions of cold 
postshock gas seen downstream in the trailing edge of the 
WCR in Fig. H», n c t ~ lO 11 ' 5 " 12 cm~ 3 s). While the ion- 
ization will not be in equilibrium in such regions, their den- 
sities are sufficiently low that their emission will be a very- 
minor contribution to the total from the system, so, for in- 
stance, it is likely that there will be no discernible effect on 
the global X-ray emission (though, as already noted, these 
regions are also the hottest parts of the WCR and thus pro- 
duce the hardest emission) . 

Figs. [S] and [S] show various slices through the simulation 
of model cwbl, which allow one to more fully appreciate the 
complex 3D structure of the winds and their collision. In 
systems with equal wind momenta, like model cwbl, the 
WCR twists like a helix out of the orbital plane. The WCR 
thickens above and below the orbital plane, and its orien- 
tation increasingly lags that within the orbital plane as one 
moves further from the orbital plane. The dense gas remains 
clumpy and fragmented, and is gradually accelerated as it 
moves out of the orbital plane. 

Fig. [7J a ) shows histograms of the mass within 120 R© 
of the system centre of mass, divided into density bins of 
width 0.1 in log. Gas at all temperatures is included, but the 
densest gas within 1.2 stellar radii is excluded. The highest 
density then obtained is 2.4 x 10 _12 gcm -3 , which occurs 
in the cold postshock gas between the stars. The mass as a 
function of density rises rapidly as the density decreases until 
p w 2 x 10~ 14 g cm~ 3 , after which it is roughly constant until 
the density is approximately 2.5 dex lower. A low density tail 
extends to p ~ 2 x 10~ 18 g cm -3 . Such rarefied gas exists in 
downstream regions of the WCR which are "shadowed" by 



the cool dense ribbon of gas closer to the WCR apex (cf. 

Fig. m . 

To better understand Fig. [TJa) we also plot the mass 
distribution between radii of 12 — 120 R© from a single 06 
star. If material leaves the surface of the star at terminal ve- 
locity the slope of the resulting mass distribution (labelled as 
"vinf" in Fig.[3i) is —0.5. An accelerating wind induces cur- 
vature into this slope, such that there is relatively more mass 
at higher densities (this distribution is labelled as "CAK" 
in Fig. [7£i). Comparing the "CAK" and "cwbl" profiles we 
see that the overall effect of a WCR is to create much more 
higher density gas in the immediate circumstellar environ- 
ment than would otherwise be the case. 

Fig. Efb) shows the mass within radii 12 — 120 Rq 
as a function of speed, in bins of 50 kms" 1 width. The 
"CAK" profile shows that most of the mass in the single- 
star case has been accelerated to nearly the terminal wind 
speed. In contrast, in model cwbl the distribution is rea- 
sonably flat between 1000 — 3000 km s -1 , with less mass at 
lower speeds. Fig. [T^b) also reinforces the point that there 
is significant mass moving at speeds up to 500 kms -1 faster 
than the terminal speeds of the individual stars, while a 
smaller proportion of mass exceeds the terminal wind speed 
by ~ 1000 kms" 1 . 

Fig. [Tfc) shows the mass in model cwbl as a function 
of temperature, in bins of width 0.1 in log. Most of the mass 
(unshocked plus shocked material) is cold (f» 10 4 K). Some 
material is shock heated above 10 7 K, with the majority 
of the hot (T > 10 5 K) gas at temperatures of a few mil- 
lion degrees. Note that the mass distribution between a few 
times 10 4 K and a few times 10 7 K bears some resemblance 



10 J. M. Pittard 




500 1000 1500 2000 2500 3000 3500 4000 
Speed (km s" 1 ) 




10° 

Temperature (K) 



cwb2 
cwb3 
cwb3o6 
cwb3o8 




500 1000 1500 2000 2500 3000 3500 4000 
Speed (km s" 1 ) 




10° 

Temperature (K) 




10° 

Temperature (K) 



Figure 7. Histograms of various quantities in the models. The left column shows data from models cwbl and cwb2 within a radius 
of 120 R,0 of the system centre of mass, while the right column shows data from models cwb2 and cwb3 within 285 Rq of the centre 
of mass. Plot a) shows mass versus density, plot b) mass versus speed, plot c) mass versus temperature, and plot d) volume versus 
temperature. Plot e)-h) show the equivalent analysis of models cwb2 and cwb5. 



to an inverted "cooling curve" , with less gas at temperatures in Fig. Eld). Most of the volume within the investigated 
where cooling is very rapid. This behaviour is quite different region contains cold gas. 
to the ma ss distributions of adi abatic wind-wind collisions 
plotted in lLemaster et al.l (|2007r ). 

The volume of gas as a function of temperature is shown 
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Figure 8. Density (top), temperature (middle), and speed (bottom) slices of model cwb2 in the orbital plane (left), and below the orbital 
plane (middle and right). The colour scales are the same as in Fig. [5] 



3.2 Model cwb2 

A density snapshot of model cwb2 is shown in Fig. QJb) 
where it can be seen that the character of the WCR is com- 
pletely different to that of model cwbl. The most strik- 
ing difference is that the postshock gas remains largely 
adiabatic, flowing out of the system while still relatively 
hot. This was anticipated since the estimated value of \ 
is substantially larger than unity. The greater distance be- 
tween the stars allows the winds to accelerate along the 
line-of-centres to ~ 1630 km s" 1 before their collision (c.f. 
a speed of 1820 km s" 1 at the same distance in the sin- 
gle star case). This results in postshock temperatures up 
to « 4 x 10 7 K (Fig. Eg), while the maximum postshock 
density is rss 4 x 10~ 15 gcm -3 (Fig. [7^), compared to 
2.4 x 10~ 12 g cm~ 3 in the cool, dense postshock gas in model 
cwbl. 

The predicted aberration angle is ~ 7°. This is again 



higher than is observed (~ 3 — 4°) for the same reasons as 
mentioned for model cwbl. Fig. QJb) also shows that there 
is a density gradient across the flow direction in the down- 
stream parts of the WCR (i.e. there is a reduction in density 
from the contact discontinuity towards the leading shock, 
and an increase towards the trailing shock). This is similar 
to the behaviour already noted in model cwbl, being caused 
by the slower exit out of the system of the gas in the WCR 
than in the winds, and the shadowing introduced by the cur- 
vature of the WCR which results in low densities ahead of 
the leading shock, and therefore seems to be a general prop- 
erty of CWBs with equal winds (the unequal winds case, 
discussed in Sec. 13.31 is more complicated). 

The typical speed of gas in the WCR as it leaves the hy- 
drodynamical grid in the orbital plane as shown in Fig. [ljb) 
is about 2500 kms -1 . In 10 days this gas travels a distance 
of 14.4 au (3100 R©). The minimum distance from the cen- 
tre of the grid to its edge is 285 R© , so we expect the WCR 
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Figure 9. Density (left), temperature (middle), and speed (right) slices of model cwb2 in the plane orthogonal to the orbital plane 
containing the stars (top), and in the plane orthogonal to this and the orbital plane (bottom). The colour scales are the same as in Fig. [5] 



to curve through 9 per cent of a full rotation, or about 33° . 
At a distance of 285 R© from the centre of mass, the WCR 
in Fig. QJb) actually displays a curvature of just over 60° , 
again indicating that the speed of gas near the apex of the 
WCR is much slower than further downstream. Hence the 
shape of the WCR near its apex is not very well approxi- 
mated by an Archimedean spiral since the speed within the 
WCR is not constant. 

We again ove r plot r esults from the dynamical model of 
IParkin fc Pittardl (|2008f ). The contact discontinuity traced 
by the red line was calculated assuming terminal wind 
speeds of 2500 kms -1 , while the green line is from a 
model with wind speeds of 1650 kms -1 . Both models as- 
sumed that the shocked wind flowed ballistically once it had 
reached 85 per cent of its terminal speed. Clearly the model 
which assumes the higher terminal wind speeds is a better 
match to the hydrodynamical calculation, which indicates 
that the downstream morphology of the WCR is not solely 
dependent on the pre-shock wind speeds along the line of 
centres. 

iLemaster et al.l (|2007l ) noted that due to the curvature 
of the WCR by the orbital motion of the stars, KH instabil- 
ities occurred even when the winds had equal speeds. There 
is no sign of this in our model, but this could be because 
of the higher ratios of Vorb/fw in some of Lemaster et al.'s 
models. The fastest wind speed attained in model cwb2 in 
the orbital plane is ~ 3200 km s -1 , near the trailing shocks 
of the WCR as in model cwbl. The lowest wind speeds in 
the orbital plane occur in the top left and bottom right of 
the panel shown in Fig. [3jb), where the speed is limited to 
~ 2400 km s -1 . Both these values are lower than the corre- 



sponding values obtained from model cwbl, indicating that 
the boosts to the acceleration of certain regions of the winds 
are greater the closer the two stars (and the centres of their 
radiation fields) are to each other. 

Fig. [4jb) displays the temperature of the WCR in the 
orbital plane. The gas near the trailing shocks cools to 
~ 5 x 10 6 K as it reaches the edge of the grid, but close to the 
leading shocks the postshock temperature is w 7x 10 5 K near 
the edge of the grid. The lower temperatures near the leading 
shocks arise because the leading shocks are more oblique at a 
given downstream distance than the trailing shock. Hence, 
in systems with identical winds the gas near the leading 
shocks is both less dense and cooler than gas near the trail- 
ing shocks, and vice- versa. The situation changes slightly in 
systems with unequal winds, as noted in Section [3.31 

In model cwb2 the postshock ion and electron temper- 
atures at the apex of the WCR equilibrate in about 2 hours. 
This is about an order of magnitude faster than the flow time 
of this gas out of the system (£fl ow ~ (D sep /2)/(v w /i) ~ 
7 x 10 4 s). Further downstream there is a thin postshock 
region where there is a noticeable delay for the electrons 
to heat to the same temperature as the ions. However, the 
thickness of these regions is typically only a few per cent 
of the width of the WCR at these positions. Given this 
and the reduced plasma emissivity compared to the apex 
of the WCR, it is not expected that there will be any ob- 
servable effects revealing T c < Ti on in systems similar to this 
model. However, slow electron heating is certainly important 
in wider syst ems, such as WR 140, as other work has alrea dy 
pointed out (|Zhekov fc Skinner|[2000l ; IPollock et al.ll2005h . 

The ionization age, n a t, is w 2 x 10 1 
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postshock gas near the apex of the WCR, and reaches 
4 x 10 13 cm -3 s at the contact discontinuity in the densest 
part of the WCR (i.e. near the stagnation point). Further 
downstream the postshock flow does not quickly reach ion- 
ization equilibrium. However, the central part of the WCR 
(close to the contact discontinuity) remains in collisional ion- 
ization equilibrium, as this gas originates from close to the 
WCR apex. 

Figs. [8] and [9] show various slices through the simulation 
of model cwb2, again to allow one to more fully appreciate 
the 3D structure of the winds and their collision. Compared 
to model cwbl, we once again see that the interaction region 
is much smoother and more stable. The top right panel of 
Fig. [5] reveals that the winds accelerate faster out of the 
orbital plane than regions within it where one of the stars 
shadows the other, due to the more effective combination of 
the stellar radiation fields. 

Fig. [T^a) shows a histogram of mass versus density 
within 120 R0 of the centre of mass of model cwb2 (exclud- 
ing material within 14 Rq of the centres of the stars). The 
distribution outlined in blue shows that most of the mass is 
at a density of ~ 5 x 10" 16 gem" 3 , and extends to ~ 4 x 
10" 17 gem -3 at the low density end and w 6 x 10 - 14 gem -3 
at the high density end, the latter representing unshocked 
wind material close to the stars. Compared to model cwbl 
the distribution is not as flat or as wide. 

Fig. Efa) also shows the mass-density distribution for 
only the hot (T > 10 K) gas. The maximum density is now 
~ 4 x 10 -15 gem -3 . The full (blue) distribution is a rough 
combination of two "CAK" distributions (one for each star) 
plus the hot (pink) distribution. However, it is clear that 
the sum of these distributions would lead to more mass at 
densities below ~ 2 x 10 -16 gem" 3 than is actually seen in 
model cwb2. The excess mass at higher densities in model 
cwb2 is because the WCR gas flowing through the spherical 
surface of radius 120 R© has p > 2 x 10" 16 gem -3 . 

The mass-speed distribution is shown in Fig.[7£b). Com- 
pared to model cwbl we find a reduced maximum speed 
attained by the gas (with the distribution turning over at 
« Uoo, rather than at « 500 kms" 1 above the single star ter- 
minal wind speed). There is also less gas moving at speeds 
below 1250 kms - , and more mass moving at speeds be- 
tween 1250 — 2600 kms" 1 . Limiting our analysis to only 
hot (T > 10 5 K) gas we see that there is a peak in the 
distribution near 1500 kms" 1 . The maximum speed of hot 
gas in the orbital plane is 2500 kms" 1 , increasing to nearly 
3000 kms" 1 out of the orbital plane. 

Fig. E^c) shows the distribution of mass versus temper- 
ature. The majority of the gas mass is predominantly cold 
(since the gas in the WCR stays hot in model cwb2 as it 
flows out of the system, the cold material represents the 
unshocked winds). However, there is also significant mass 
at temperatures ~ 10 7 K. The amount of mass at interme- 
diate temperatures (~ 10 4 — 7 x 10 5 K) should be treated 
cautiously, since it mostly represents gas in cells at the unre- 
solved shocks, and will therefore be sensitive to the numer- 
ical scheme. The bi-modal nature of the mass-temperature 
distribution of model cwb2 is also reflected in Fig. [7jd), 
which shows the volume-temperature distribution. 76 per 
cent of the volume within 120 R© of the system centre of 
mass contains the unshocked winds, while 20 per cent con- 
tains gas at temperatures exceeding 5 x 10 6 K. 



Figs. Efe)-(h) show mass and volume distributions of 
the hot gas in model cwb2 out to a radius of 285 R© from 
the system centre of mass. The distributions are similar to 
those in Figs. [7ja)-(d), but extend to lower densities and 
temperatures, and greater speeds, as expected. 

3.3 Model cwb3 

A distinct difference between the WCR in model cwb3 com- 
pared to the other models is that its apex occurs closer to 
one star (the O8V secondary star which has the weaker 
wind) than the other star (the O6V primary star) due to 
the unequal wind momentum ratio (see Fig. [T};). The gas 
downstream of the WCR apex is also pushed closer to the 
secondary star, due to the greater momentum flux of the 
primary wind. These effects, together with the fact that the 
winds are still accelerating when they collide, plus the en- 
hanced terminal wind speed of the 06V star, lead to higher 
preshock and postshock speeds in the primary wind relative 
to the secondary wind. This has three major consequences: 
i) there is a velocity shear along the contact discontinuity 
between the winds, which excites Kelvin-Helmholtz insta- 
bilities; ii) to maintain pressure balance, the preshock and 
postshock density of the secondary wind exceeds that of 
the primary wind; iii) higher postshock densities, together 
with reduced postshock speeds, results in the gas on the sec- 
ondary side of the WCR radiating more efficiently (this is 
discussed in more detail in a forthcoming paper on the X-ray 
emission) . 

There is also a clear signature of radiative inhibition: at 
the primary wind shock on the line-of-centers, the pre-shock 
speed is 1800 kms" 1 , whereas without the companion star 
it is 2010 kms - at this distance. Hence the radiation field 
of the secondary star reduces the net acceleration of the 
primary wind towards it by ~ 10 per cent. Likewise, the pre- 
shock speed of the secondary wind on the line-of-centers is 
1270 kms" 1 , whereas the single star solution has a speed of 
1450 kms -1 at this distance, indicating a reduction in the 
acceleration of « 12 per cent. Radiative inhibition thus has 
a greater effect on the seconda ry wind, as expecte d. There 
is no sign of radiative braking (jCavlev et al.lll997h . though 
this is not surprising given the low wind momentum ratio 
(hence wind opacity ratio) in our model. 

Also overplotte d on Fig. [3fc) is the comp arison with the 
dynamical model of IParkin fc Pittardl (120081 ). The match is 
again very good, with the main differences between the mod- 
els occuring in the leading arm where the contact discon- 
tinuty in the hydrodynamical model shows slightly greater 
curvature. 

Fig. [3jc) shows the speed of gas in the orbital plane. 
The maximum speeds reached in the 06 and 08 winds are 
2910 kms -1 and 2840 kms" 1 , respectively. These val- 
ues are attained near the trailing edges of the WCR at 
the boundaries of our hydrodynamic grid. More typically, 
the 06 and 08 winds reach speeds of w 2550 kms -1 and 
2300 kms -1 (a gain in directions towards the top left and 
bottom right of the hydrodynamic grid), respectively. Gas 
in the denser parts of the WCR typically reaches speeds of 
1600 kms" 1 in the orbital plane by the time it leaves the 
hydrodynamic grid. 

Using the measured pre-shock velocities along the line 
of centres to determine the cooling parameter \, we find 
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that XX ~ 28 and X2 ~ 14. Thus the postshock winds re- 
main largely adiabatic, as is indeed seen. The postshock 
winds reach maximum temperatures of ~ 5 x 10 7 K (for 
the shocked 06 wind) and ~ 3 x 10 7 K (for the shocked 
08 wind), in good agreement with theoretical expectations. 
The maximum densities attained are ~ 2 x 10~ lo gcm -3 
and 5 x 10~ 15 gcm~ 3 for the shocked 06 and 08 winds, re- 
spectively. These maximum values are all attained near the 
stagnation point. 

The ratio of post-shock densities across the contact dis- 
continuity (CD) is found to vary with downstream distance. 
The factor of 2.5 difference at the stagnation point increases 
to a factor of 7 — 10 in the arm where the 08 gas forms the 
trailing edge, but becomes a factor of 3 or so in favour of 
the 06 wind (i.e. the shocked 06 material is denser than the 
shocked 08 material) in the opposite arm. The temperature 
difference across the CD increases in the WCR arm where 
the 06 material is at the trailing edge, whereas it inverts 
(i.e. the shocked 08 material becomes hotter, through less 
rapid cooling, than the 06 material) in the opposite arm. 

The postshock electron and ion temperatures rapidly 
attain equilibrium in all parts of the flow, in common with 
the previous models. The significantly higher density and 
lower temperature of the shocked 08 wind at the stagnation 
point leads to more rapid ionization equilibrium compared to 
the shocked 06 wind. At the WCR apex, the ionization age 
of the shocked 06 wind is typically 4.5 x 10 12 cm -3 s, while 
nearer to the CD it increases to « 3.5 x 10 13 cm~ 3 s. The 
ionization age of gas near the CD further increases down- 
stream, as already noted for the other models. 

Fig. [7J e ) shows the mass-density distribution of hot gas 
(T > 10 5 K) in model cwb3 within a spherical volume cen- 
tered on the system centre of mass and extending to a ra- 
dius of 285 R© (again excluding material within 14 R© of 
the centres of the stars). The distribution is similar to that 
from model cwb2 over an identical volume, though the total 
mass is somewhat less owing to the smaller mass-loss rate of 
the secondary star in model cwb5. However, it is interesting 
to find that the mass at higher densities is dominated by 
material from the 08 wind, and vice-versa. This is because 
the shocked 08 wind is denser at the WCR apex than the 
06 wind (due to the latter's greater pre-shock speed), and 
remains denser as it flows out of the system due to the cur- 
vature/confinement of the 08 gas as the WCR bends behind 
the 08 star. 

The mass-speed distribution is shown in Fig.[?Kf). Com- 
pared to the distribution obtained from model cwb2, the 
cwb3 distribution extends to a slower maximum speed. This 
is due to the lower terminal velocity of the 08 wind, and 
the reduction in the combined radiative driving force due to 
the reduced luminosity of the 08 star. Model cwb3 has more 
mass at speeds < 1400 km s -1 - this is dominated by mate- 
rial from the 08 star, and indicates the slower acceleration 
of its wind, and subsequently the lower speed of its post- 
shock gas. Conversely, material from the 06 star dominates 
the hot gas moving at speeds above 1600 km s -1 . These dif- 
ferences may have implications for the emission line profiles 
from the WCR of such systems. 

Fig. [TJg) shows the mass-temperature distribution. 
Compared to model cwb2, the distribution from model cwb3 
contains less mass overall (as already noted), but also has 
gas at a higher maximum temperature, since the 06 wind 



has more room to accelerate along the line-of-centres before 
being shocked. Model cwb3 also contains much more mass 
at T < 2 x 10 6 K because of the slower acceleration and 
speed of the 08 wind. 

The volume-temperature distribution of hot gas (T > 
10 s K) is shown in Fig. [TJh). Compared to model cwb2 we 
find that the volume of hot gas above 4 x 10 7 K is greater 
in model cwb3, due to the faster pre-shock speed of the 
06 wind. The volume of hot gas is slightly less between 
4 x 10 6 - 4 x 10 7 K, is comparable at rts 3 x 10 6 K, and is 
greater for T < 1.4 x 10 6 K. Shocked 06 gas dominates the 
volume of hot plasma at T > 2 x 10 6 K, below which the 08 
wind dominates. 

The emission measure (n 2 V) of the hot gas, divided 
into temperature bins of width 0.1 in log, peaks at 9.2 x 
10 54 cm' 3 at T « 1.75 x 10 7 K (this is not shown in 
Fig. UJ. The shocked 08 wind is almost completely dom- 
inant, with the emission measure from 06 wind material 
peaking at 2.4 x 10 54 cm -3 at the slightly hotter tempera- 
ture of T w 2.25 x 10 7 K. Thus the shocked 08 wind should 
dominate, for example, the thermal radio and X-ray emis- 
sion generated in the WCR - in fact, we find that the 08 
wind dominat es the intrinsic X-ray emission below 4.5 keV 
(Pittard 2009, in preparation). The total emission measure 
is 7.1 x 10 55 cm -3 , with the shocked 08 wind contributing 
70 per cent of this value. Since the shocked 08 material in 
the leading arm of the WCR is compressed into about half of 
the volume of the shocked 06 material in this arm (note the 
high density contrast in Fig. []J), it is interesting to exam- 
ine the contributions of each wind to the emission measure 
from each arm of the WCR. We find that the shocked 08 
material dominates the emission measure in the leading arm 
(2.8 x 10 55 cm~ 3 versus 0.8 x 10 55 cm~ 3 for the shocked 06 
wind) . In the trailing arm the ratio of the emission measures 
between the winds is much closer to unity (2.0 x 10 55 cm~ 3 
for the shocked 08 wind, versus 1.4 x 10 55 cm -3 for the 
shocked 06 wind) , though the 08 material again dominates 
the emission. 

3.4 Model cwb4 

As one might expect given the large differences between 
models cwbl and cwb2, the structure and properties of the 
WCR in model cwb4 undergoes rich and dramatic changes 
as a result of the varying separation between the stars as 
they progress in their respective orbits. Fig. [TU] shows the 
evolution of the mass density in the plane of the orbit. At 
periastron the WCR is highly radiative and subject to strong 
dynamical instabilities. The cooled postshock gas is sheet- 
like, as in model cwbl, and maintains a coheasive structure. 
As the stars separate and the cooling of the postshock gas 
becomes less severe, this sheet fragments and discombobu- 
lates into smaller, individual clumps. By apastron, the WCR 
is largely adiabatic. At both periastron, when the stars are 
separated by the same distance as in model cwbl, and at 
apastron, when the stars are separated by the same distance 
as in model cwb2, the wind-wind interaction somewhat re- 
sembles the WCR in these other models (cf also Fig. I15[) . 
For instance, the velocity profiles between the stars along 
the line of centres from model cwb4 at periastron and apas- 
tron are in good agreement with those from models cwbl 
and cwb2, respectively. 
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Figure 10. Density snapshots of model cwb4 at specific phases in the orbit, which progresses from left to right, and from top to bottom. 
The top left panel shows the model at periastron, and subsequent panels increase the orbital phase in steps of 0.1. The right most panel 
of the 2nd row from the top shows the stars at apastron, with periasron again occuring in the middle panel of the bottom row. The 
colour scale is logarithmic, spanning 10 -19 gem -3 (blue) to 10 — 11 gem -3 (red). The panels have sides of length 240 Rq. 



Nevertheless, there are some interesting differences. For 
instance, some dense clumps of cool gas exist within the 
WCR even at apastron, in contrast to the entirely hot state 
of the WCR in model cwb2. These clumps were formed dur- 
ing the previous periastron passage. Because of their high 



density and inertia relative to the surrounding hot, more 
rarefied, plasma, it can take some considerable time for the 
faster flowing hot plasma to push the clumps out of the sys- 
tem. Clearly this timescale is of the order of half the orbital 
period of the stars in model cwb4. The slow acceleration 
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Figure 11. As Fig, llOl but showing temperature snapshots. The colour scale is logarithmic, spanning ^ 10 5 K (white) to 10 s K (pink). 



of the clumps out of the system also results in some of the 
clumps crossing the global shocks bounding the WCR and 
moving into the fierce high speed environment of the un- 
shocked winds (this is clearly seen at phases 0.4 — 0.6 in 
Fig. I10[) . This process occurs because the hot plasma in the 
WCR responds quite quickly to the movement of the stars 
(and thus rotates with the stars), whereas the direction vec- 
tors of the dense clumps respond much less rapidly and the 



clumps tend to maintain straighter paths. For this reason the 
clumps exit only out of the trailing shocks of each arm of the 
WCR, and not out of the leading shock. Once formed, the 
clumps are gradually destroyed by thermal evaporation and 
ablation, with larger and denser clumps surviving longer. 

The presence of clumps in or around an otherwise adia- 
batic WCR may affect some observational signatures. For in- 
stance, gas mixed from the clumps into the surrounding flow 
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Figure 12. Density (top), temperature (middle), and speed (bottom) slices of model cwb4 parallel to the surface of the WCR near its 
apex at periastron. The left panels pass through the WCR downstream of its apex, the middle panels pass through the centre of one of 
the stars, and the right panels pass through the stagnation point and apex of the WCR. The colour scales are the same as in Fig. [5] 



"mass- loads" it, making it denser. The effect of the clump on 
the temperature of the surrounding flow is more complex. A 
bowshock forms upstream when the surrounding flow is su- 
personic with respect to the clump, heating the surrounding 
flow. This is seen, for example, when the clumps move out 
of the WCR and into the unshocked winds. Bowshocks also 
form upstream of clumps within the WCR when the clumps 
are far enough downstream that the hotter flow past them is 
supersonic - this is the case for the lowermost clump above 
the WCR apex in the apastron panel of Fig. 1101 Friction be- 
tween the clump and a surrounding subsonic flow also heats 
the flow. Further downstream, the addition of mass into a 
flow generally leads to lower temperat ures once the injecte d 
mass is fully mixed into the flow (e.g. IPittard et al, 2003a). 
though this will require reconnection of field lines in a mag- 
netic environment. The presence of clumps outside the main 
WCR may also lead to more of the winds' kinetic power be- 



ing processed through shocks than would otherwise be the 
case, and may affect the variation of, for example, the X-ray 
emission with phase. Finally, the boundary layer between 
clumps and the surrounding diffuse flow contains intermedi- 
ate densities, temperat ures, and velocities, all of which affect 
the resulting emission (jHartauist fc Dvsonlll993h . Line emis- 
sion is expected to be particularly sensitive to the conditions 
in this region, and such regions may help to expla in dis- 
crepancies in current models (e.g. Irlenlev et al.ll2008l ) where 
mixing of cold material into a hotter medium is neglected. 

Another difference with our previous models is that cold 
dense gas does not exist far downstream of the apex of the 
WCR when the stars are at periastron, whereas it does in 
model cwbl (compare the top left panels in Figs. HOI and llll 
with Figs. [Ha) and|4ja), and remember that the stellar sep- 
arations are identical). This is simply because the down- 
stream gas in model cwb4 at this phase reflects the less 
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Figure 13. 3D temperature surface from model cwb4 at apas- 
tron. The surface traces gas at a temperature of 3 X 10 6 K. 

radiative environment which existed within the WCR when 
the stars were more widely separated. 

Not surprisingly, there is a large variation in the aber- 
ration and curvature of the WCR in model cwb4 with or- 
bital phase. The aberration angle reaches a maximum at 
periastron due to increased stellar velocities and reduced 
preshock wind speeds. At periastron, v OT b = 334 kms -1 , 
while u w w 710 kms , giving an expected aberration an- 
gle of « 25°, in good agreement with the measured angle 
of ~ 21°. At apastron, v OI b drops to 156 kms -1 , while v w 
increases to ~ 1665 kms -1 , giving an expected aberration 
angle of m 5°, again in good agreement with the measured 
angle of ~ 4°. Note that the orbital speeds at periastron 
and apastron are different to the speeds in models cwbl and 
cwb2, so the aberration and curvature of the WCR differs 
compared to the circular orbit models (see also Table [IJ. 

Fig. [11] shows the temperature in the orbital plane of 
model cwb4. The hottest gas occurs at apastron (T max m 
5 x 10 7 K), and its temperature is similar to the hottest gas 
in model cwb2. While the maximum gas temperature attains 
a minimum around periastron (T max ~ 2 x 10 7 K), this value 
is found in localized bowshocks enveloping dense clumps in 
the downstream flow and is sensitive to rapid changes in the 
local dynamics. It also relates to a very small fraction of the 
total mass within the WCR, and hence is not representa- 
tive of the typical temperatures which are much lower (cf. 
Fig. 115b ). Fig. 1111 also reveals that at periaston the trailing 
edge of each arm of the WCR is hotter than the leading 
edge. This is because of the severe curvature of the WCR 
due to orbital motion, which results in the shock at the trail- 
ing edge staying more normal to the preshock wind vector. 
This effect is not seen at apastron when the effects of or- 
bital motion are less intense (see Fig. Hip , However, model 
cwb2 suggests that gas at the trailing edge of each arm of 
the WCR will become hotter than gas at the leading edge at 
distances further downstream than can be probed in model 
cwb4 (model cwb2 has a larger grid - see Table [3]) . 

Fig. [12] shows various slices through the 3D hydrodyan- 



mical grid when the stars are at periastron, again to allow 
appreciation of the 3D structure of the winds and their col- 
lision. The slices in the right-most panels pass through the 
stagnation point at an angle aligned with the shock surfaces 
between the stars, and so sample a relatively large distance 
through the WCR despite its small width. This figure il- 
lustrates again that cold dense gas exists near the apex of 
the WCR at periastron, but not further downstream at this 
phase. 

Fig. [13] displays a 3D temperature surface from model 
cwb4 at apastron. The sheet-like structures of the global 
shocks bounding the WCR are clearly visible. However, this 
figure also highlights that while cold gas has been largely 
cleared out of the central regions of the WCR, some still 
remain further downstream. The large bowshock projected 
towards the bottom right of the plot envelopes the dense 
blob of gas seen outside of the WCR towards the bottom 
left of the (j> = 0.5 panel in Fig. [TO] Fig. [T3] shows tem- 
perature surfaces viewed from above the orbital plane at 
phase (j> — 0.4. Fig. I14f a) traces gas at a temperature of 
2 x 10 J K, which highlights the global shocks bounding the 
WCR and the interface between cold, dense clumps embed- 
ded within the hotter, more rarefied plasma of the WCR. 
Figs. I14f b) and (c) show surfaces of progressively higher 
temperatures. As the temperature of these surfaces increases 
the intermediate-temperature material ablated from the cold 
clumps is revealed. Comparison between Figs. [TjJ] and I14f c'l 
highlights the gradual removal of cold dense clumps from 
the WCR as the orbit progresses away from the previous 
periastron passage. 

Fig. I15f a) shows the mass within 120 Kq of the sys- 
tem centre of mass (excluding the stars) as a function of 
density from model cwb4 at periastron, and from model 
cwbl. The distributions are reasonably similar, with the 
main differences being a relative excess of mass at high 
(p > 5 x 10 -14 gem -3 ) and intermediate (p ~ 10~ 15 gcm~ 3 ) 
densities, and a lack of low density material (p < 3 x 
10 -16 gem -3 ) in model cwb4. The first of these is explained 
by the recent history of the interaction in model cwb4 - i.e. 
that there is hotter gas in the downstream regions of the 
WCR because this gas was shocked and advected out of 
the central regions of the system before cooling became as 
strong as it is now that the stars are at periastron. The last 
is explained by a reduction in the degree of "shadowing", 
since prior to periastron the rate at which the stars sweep 
out angle is much smaller than the rate in model cwbl. 

A comparison between the mass versus density his- 
tograms of models cwb2 and cwb4 at apastron (Fig. 115b ) 
again reveals similar distributions, with the most notice- 
able difference being a relative excess of mass at p > 5 x 
10 -lo gcm -3 in model cwb4. This now reflects the presence 
of cold dense clumps in the WCR from the previous perias- 
tron passage. A comparison of the mass of hot (T > 10 5 K) 
gas reveals a similar excess at higher densities, due to the 
intermediate temperatures of gas ablated off the cold dense 
clumps. 

The mass versus speed histogram of models cwbl and 
cwb4 at periastron are shown in Fig. I15f c). The material 
in model cwb4 does not quite reach the maximum speeds 
seen in model cwbl, most likely because the combined ra- 
diative driving force, which is greater the closer the stars 
are to each other, would have been weaker in model cwb4 
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prior to periastron. Model cwb4 also displays an excess 
of mass relative to model cwbl at speeds < 1400 km s" 1 . 
Fig. I15f d) compares the distributions between models cwb2 
and cwb4 at apastron. They are remarkably similar, except 
again model cwb4 displays a relative mass excess at speeds 
< 1400kms _1 , which is likely due to the presence of cold 
clumps. The distributions of hot (T > 10 5 K) gas are also 
remarkably similar. 

Fig. I15f e) shows mass versus temperature histograms 
from models cwbl and cwb4 at periastron. Model cwb4 has 
more than twice as much mass at temperatures above 10 6 K, 
but less mass at lower temperatures, due to the reduced 
cooling in the post-shock history of the downstream gas in 
model cwb4. At apastron, model cwb4 has slightly less mass 
at T > 8 x 10 6 K compared to model cwb2 (Fig. [151), but 
significantly more mass between temperatures of 4 — 8 x 
10 6 K. 

The volume versus temperature histograms from mod- 
els cwbl and cwb4 at periastron are shown in Fig. I15f g). 
There is a significantly greater volume of gas at T ~ 10 6 K 
in model cwb4. This excess combines with a reduction in 
the relative volume of gas at intermediate temperatures 
(1.5 x 10 4 < T < 6 x 10 5 K). These differences reflect the 
hotter thermal history of the shocked gas in model cwb4 
prior to periastron. At apastron in model cwb4, the vol- 
ume of hot gas at temperatures above 8 x 10 6 K is 25 — 60 
per cent smaller than in model cwb2. Model cwb4 generally 
also has smaller volumes of gas at intermediate tempera- 
tures (1.5 x 10 4 < T < 4 x 10 6 K) - this is likely a numerical 
effect due to the higher resolution in model cwb4 compared 
to model cwb2, which leads to a narrower layer of intermedi- 
ate temperature cells at the unresolved shocks. Model cwb4 
has almost an order of magnitude greater volume of gas at 
5 — 6 x 10 6 K than model cwb2. Such temperatures are en- 
countered in the mixing regions downstream of the clumps 
in model cwb4 (c.f. Fig. [14)) . 



4 SUMMARY AND CONCLUSIONS 

The model results described above reveal a wealth of in- 
teresting dynamical effects which occur in short period, 
O+O binaries, including radiative inhibition, rapid post- 
shock cooling, and powerful instabilities, all of which had 
been expected. However, there have also been some sur- 
prising findings, such as the ability of cold dense clumps 
to exist for substantial periods of time in the WCR when 
current postshock conditions are largely adiabatic. Another 
surprise is that the slow exit velocity of these clumps out of 
the system allows some of them to exit the WCR through 
its trailing shock. More generally, the presence of clumps 
increases the complexity of the WCR, affecting the density, 
temperature and velocity of the shocked plasma, and al- 
lowing additional material to be processed through normal 
shocks. As a result, the emission from the WCR may be 
affected in a non-trivial manner. 

In future work we will study the multi-wavelength emis- 
sion from these models, and will repeat our study for WR+O 
systems where we expect radiative braking to have a signif- 
icant effect on the dynamics, and the observed emission to 
show greater orbital modulation due to the larger differences 
in the attenuation through the stellar winds. Braking should 



also be important in O+O systems where there is a larger 
ratio of wind opacities than explored here. We will also turn 
our attention to modelling specific systems, with the aims of 
determining key system parameters and understanding the 
natures of the diverse range of wind- wind interactions which 
occur. 
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Figure 15. Histograms of various quantities in model cwb4, compared to models cwbl and cwb2. The left column shows data from 
models cwbl and cwb4 with the later at periastron, while the right column compares model cwb2 with cwb4 at apastron. All histograms 
examine the material within a radius of 120 Rq of the system centre of mass, and exterior to 1.2 stellar radii around the stars. 



